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An equal time version of odd-frequency pairing for a generalized t — J model is introduced. 
It is shown that the composite operators describing binding of Cooper pairs with magnetization 
fluctuations naturally appear in this approach. The pairing correlations in both BCS and odd- 
frequency channels are investigated exactly in ID systems with up to 16 sites. Our results indicate 
that at some range of parameters odd-frequency correlations become comparable, however smaller 
than BCS pairing correlations. It is speculated that the spin and density fluctuations in the frustrated 
model lead to the enhancement of the odd gap susceptibilities. 
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The odd- frequency gap (odd gap hereafter) originally 
proposed for spin triplet pairing by Berezinskii |l|, and 
recently extended on a spin singlet pairing ||, was the 
subject of a recent investigation using the Monte Carlo 
method on the Hubbard model ||. The motivation for 
the current study was the observation that odd pairing 
can be a possible pairing channel in strongly correlated 
systems M. 

As it was shown in (5) , the odd gap can be represented 
in another form as an equal-time (i.e. even in frequency) 
composite four fermion operator constructed out of equal 
time pairing field operator and a magnetization operator 
(generally a particle- hole operator). The total charge 
of this four fermion composite operator, i.e. the charge 
with which this object transforms under electromagnetic 
U(X), remains q = 2e as it should. This representation 
offers further insight into the nature of condensate of odd 
pairing. 

To explore the possible existence of the odd-frequency 
gap SC condensate in strongly correlated systems, we 
compared quantitatively the odd and BCS pairing cor- 
relations in the recently proposed ID t — J model in a 
staggered h field . 

One of the important advantages of studying ID vs 
2D models is that in ID case we can numerically in- 
vestigate much larger systems and thus get information 
about pairing correlations at larger distances. For exam- 
ple, in this study we report on a numerical calculation 
of both BCS and odd pairing correlations in ID system 
with up to 16 sites. This should be compared with a 
typical Monte Carlo calculations in 2D of 6 x 6 sites. 
Large system size allows us to find pairing correlations 
at much greater distances. To make the physics of this 
particular ID model closer to real 2D physics we imposed 
the external staggered magnetic filed. It has been shown 
in that the external staggered field h in ID simu- 
lates some 2D effects otherwise not characteristic for or- 
dinary ID t — J model. More precisely, it induces a longer 
range antiferromagnetic ordering which leads to the spin- 
string phenomenon B, strong mass renormalization Q 
and formation of bound hole pairs JTo| . The phase dia- 
gram of the t—J — h model shows three distinct regions 
@. When J < J c (h) there are no bound states, when 
J c (h) < J < J s (h) system consists of bound hole pairs 
and when J > J c (h) phase separation into a hole- rich 
and spin-rich phase occurs. Furthermore, since h ^ de- 
stroys the spin-rotation invariance of the system and in- 
duces a gap in the spectrum of magnetic excitations, the 
model belongs to the Luter Emery universality class jllj . 
The charge exponent K p , which determines the behavior 
of correlation functions at large distances, was proven to 
be K p > 1 in a wide, physically relevant regime of pa- 
rameters, indicating on dominant SC correlations 

In this paper we will report on a comparative study of 
even (i.e. BCS) and odd pairing correlations in spin sin- 
glet and triplet channels in ID t—J— h model using 
Lanczos exact diagonalization method. First, we will dc- 

rWp an prmal timp vpreinn nf nrlrl cran anH Qnnw that* 



equal time order parameter describes the condensai 
composite operator, which is the product of equal 
pair field and magnetization operator; b) the equal 
version of spin triplet odd field operator is even (P = 
under parity transformation of relative coordinates ai 
represented by a product of both spin singlet and tr: 
equal time pairing field and appropriate componen 
magnetization operator (see below) . The spin singlet 
operator can be written as a product of magnetizs 
operator and equal time spin triplet pairing fields 
is odd (P = —1) under spatial inversion. These feat 
agree with general symmetry requirements for odd sir 
and triplet gaps J|J§|. Next, we compare pairing cor 
tions in BCS and odd channel. The BCS correlations 
always large in comparison with odd correlations in 
t — J — h model at small hole-doping. However, in 
extended t — t' — J — h model with next nearest n« 
bors hopping amplitude t', the odd pairing correlat 
increase and become comparable with the BCS pai 
correlations. 

The main conclusion of this study is that in the ; 
t — J — h model the odd gap correlations are generii 
small. However, in the extended t — t' — J — h m 
these correlation are increased substantially (by on< 
der of magnitude). Moreover they become of the s 
order as the largest BCS correlations for a particular 
of t' . From our numerical calculations it also unaml 
ously follows that the odd pairing correlation lengl 
increasing rapidly with increased hole doping, in conl 
to suppression of the BCS correlation length. This c 
indicate that odd pairing becomes the dominant pai 
channel close to phase separation region, where par 
and spin density fluctuations are strongly enhanced 

We start with the hamiltonian of a generalized 
t — J — h model |(| in the following form 

H = -f^(ct C i+liQ! +ff.C.)-i'^(c] > aCi+2 ) a + J 

iot tot 

+ jJ2(Si ■ S i+1 - fam+i) -hJ2(-iySi, 

i i 

where c~i >a — Ci iQ ,(l — rii^ a ) is projected fermionic 
erator onto the space of singly occupied sites, rii is 
corresponding fermion number operator, S% is a spir. 
erator, h is the strength of the external staggered f 
and the t' term represents the next-neighbor hopping 
simplify our notation we omit in the rest of the papei 

sign above creation and annihilation operators. 

First we consider a general two-particle gap funct 

A Q/3 (k, 2t) = (r T c Qik (r)c /3i _k(-T) ± (r -> -t)}, 

where T T is Matsubara time ordering operator and w> 
sume pairing in the center of mass momentum zero, 
constructed the gap function in (||) to be explicitly > 
(+ sign) or odd (— sign) in the imaginary time r. 
recently suggested HQ] , the only requirement for pai 
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is even under simultaneous space and imaginary-time re- 
versal, i.e. k — > — k, t — > — t. Consequently, in the 
case of the singlet pairing, the gap function can be either 
even in both space and imaginary-time, leading to the 
conventional BCS gap, or can be odd in both space and 
imaginary-time, leading to the odd gap. Furthermore, 
in the case of triplet pairing the gap function must be 
odd under simultaneous space and imaginary-time rever- 
sal again leading to two solutions: a conventional BCS 
where gap function is even in r and odd in k or odd gap, 
which is odd in t and even in k. 

We are mainly interested in the lowest order contribu- 
tions to the gap in the relative pair-field time 2r. Then, 
assuming analyticity of A a p in small t, we have 

K7 n = A l°j + 0(r 2 ) (3) 
A° a f = 2rA^ + 0(T 3 ), (4) 

where the equal-time contribution in Eq.(Q) corresponds 
to the even gap, which is the conventional BCS gap, 
whereas the higher order corrections in r generate fur- 
ther multi-particle dressing. By definition there is no 
zcroth order term in r (i.e. equal-time pairing) term in 
the case of odd gap Eq.(Q). To calculate the lowest equal 
time pairing field associated with odd-pairing we take the 
time derivative of the gap function 

dA a/3 (k, 2r) _ A (i), km _ 
^ | T =0 - A Q/3 (k,U) - 

(c Q ,kC/5,_k - Ca,kC/5,-k), (5) 

where c 7jP , c 7iP are calculated at r = 0. 

The odd-r pair-field operator can be written in the 
following form 

A oc ^2 A[ i 1 j(k)(c Qjk C/5,-k - c aik c/5,_k)- (6) 

k 

We are still free to choose A^2(k) to be an even (odd) 
(P = ±1) under space transformation k — > k, thus 
obtaining a triplet (singlet) odd gap pair-field operator. 
To rewrite Eq. (^) in a more explicit form we proceed 
by calculating equations of motion for a Hamiltonian of 
t — t' — J — h model in the r-space 

Ci, a = —Mia a /3C itl 3 - h{-l) l a z aa c i>a 
+t(Ci-l,a + <H+l,a) + t'(ci- 2 ), (?) 

with Mi = J(Si—\ + Si+i) and a a p are Pauli matrices. 

Assuming that the gap function Al^ (k) is sin k in the 
case of odd singlet pairing and cos k for the S z = triplet 
and that Aj* (k) = cos k for S x = 1 triplet, we obtain af- 
ter a Fourier transformation from Eq. (]?]) and Eq. (||) the 
following compact form of odd singlet and triplet pairing 
operators: 



KtngleM) OC + S %+2 ) (<7 B *) o/J <*,<,<*+!,, 

A triplet,S z =o(n) OC (Si-1 (0- V a) ap - S l+2 {(7 V ' B) ^ 

^ P iet,s z =±i(n) oc (Si-! (a ± a z a) ap - S l+2 {3 ± < 

In deriving Eqs. (^,||,^0|) the terms proportional tc 
and h drop out, since they contribute only to the 
pairing operators. 

It is instructive to investigate the composite struc 
of odd pairing operator from Eqs. (PJ9|,^0|) more clo 
Consider A°f^ glet , which is a product of magnetize 

operator (Si— i + £^+2) with the BCS spin triplet ] 
ing field A^ et oc (cr y a) a/} c ljQ c i+ i^. Since produ< 
these two operators form a scalar in spin space, it in 
diately follows that A°f^ glet is indeed a singlet. The 
tial parity of composite operator is the product of p; 
of magnetization operator P ma gn = +1: at i — 1 <-» i 
we have (Si-i + Si+2) (<Si-i + Si+2) and the parii 
^triplet' which is P bcs = -1- Thus it follows that 
composite odd operator from Eq. (||) represents the 
parity P = — 1 spin singlet pairing filed. The analo. 
consideration can be done for odd triplet operator as ' 

From the structure of compi 
operators Eqs. (||,[|[l0|) it is clear that the odd pai 
correlations will be suppressed due to phase spaa 
strictions. However, there is an important caveat: < 
to the phase separation region in the t — J model, 
and density fluctuations are strongly enhanced at u 
mediate length scale . A snapshot of the system 
show metallic regions with a large density of holes 
antifcrromagnetic regions with less holes. These flue 
tions are soft and their spectral density A(u>) has a st: 
peak at small frequencies. Under these circumstance! 
soft spin boson attached to the Cooper pair in the ( 
posite operator is readily available in the system anc 
phase space restrictions become less important. 

From this argument we should expect that the 
pairing susceptibilities will increase with frustratio 
t — J model. We indeed find some enhancement of 
pairing correlations in a frustrated t — t' — J mode 
we will show below. 

We will now turn to our numerical results. We so 
the t — J — h model on a chain with periodic boi 
ary conditions using the standard Lanczos techni 
We restricted our calculations to the systems with . 
16 N h = 2, N = 14 N h = 4 and N = 12 N h = 6, w 

represents the number of sites and N% the numb< 
holes in the system. In all cases investigated the grc 
state has the quantum numbers S z = 0, k = 0. I* 
that due to the staggered external field the total sp 
no longer a good quantum number. 

We searched for the most favorable pairing channe 
ing the equal-time pair-field susceptibility in the folio 1 



form 



P(r) 



(ID 



where we took for A(r-j) the odd gap, Eqs. (p|,P|JlC|) or 
BCS gap pairing fields. We used a standard definition for 
the BCS pairing fields: A BCS (ri) oc Cj^cj^+i ic^iCj^+i 
where - (+) represents singlet (triplet) S z — pairing 

field and A^2,s,=i( r ») = c 1\« c T,i+i is a triplet 5 Z = 1 
pairing-field. 

In Fig. (1) we present density-density correlation func- 
tions g(r) = J2i( n h,inh,i+r), where nh,i = 1 - n*, in the 
system of N = 14 sites and 2V/, = 4 holes at /i/i = 1.0 
and two different values of J. At J/t = 0.5 (open circles) 
g(r) exhibits a peak at the largest possible distance tak- 
ing into account periodic boundary conditions r rn = N/2 
indicating that no bound hole-pairs are in the system. 
However, at a greater value of J/t = 2.25 there are two 
peaks in g(r) at r = 1 and r = r m which are consistent 
with the picture of two separate bound hole-pairs. 

In Fig. (2a,b) BCS (open symbols) and odd gap (filled 
symbols) pairing susceptibilities P are shown as functions 
of distance r between pairs for the same choice of param- 
eters as in Fig. (1). At J/t = 0.5, Fig. (2a), where there 
are no bound pairs in the system, pairing susceptibilities 
rapidly decay with the distance. This rapid decay is also 
in agreement with the value of the charge exponent K p 
being less than unity 0. In Fig. (2b) we present pairing 
susceptibilities at J/t = 2.25. We observe strong en- 
hancement of the BCS pairing susceptibilities at larger 
distances, however, there is no major enhancement in 
the odd gap channel. The enhancement of the pairing 
susceptibilities is consistent with the formation of bound 
hole-pairs, enhancement of the charge exponent K p > 1 
and the proximity of the phase separation . 

In our search for a possible enhancement of odd 
gap pairing susceptibilities we also included a second- 
neighbour hopping term t'/t. The effect of the next- 
neighbor hopping is strongly asymmetric. Negative 
t'/t < leads to: a substantial spin-charge coupling 
(nearly absent in the ID t — J model), formation of fer- 
romagnetic polarons and finally phase separation Jjl|. 
Since finite h > enhances antiferromagnetic type of 
spin ordering, we argue that t'/t < induces intability 
in the t — t' — J — h model. On the contrary, the ef- 
fect of t'/t > on the pairing correlations is much less 
pronounced. In Fig. (2) we present the comparison of 
largest pairing susceptibilities in the BCS and odd gap 
channel at J/t — 1.5, h/t — 1.5 and two different val- 
ues of t'/t = 0, —0.3. Clearly there is an opposite, weak 
effect of t' on pairing susceptibilities. On the one hand, 
the BCS pairing correlations diminish when t' is switched 
on; on the other hand, we observe small but consistent 
increase of odd gap susceptibility at large distances. 

Finally, we investigated the effect of doping on pair- 
ing in odd and BCS channel. To present results in a 
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resenting a sum over pairing off-diagonal susceptibil 
Ps = J^r l-P( r )l/-^; where we choose r > 1 to a 
overlap of nonlocal operators with themselves. ( 
sider also the first moment of P(r) which is del 
as £ = J2 r r\P(r)\/N J2 r \P(. r )\ aim gives the cor 
tion length of superconducting fluctuations. In T 
we present results for Ps and £ for three different 
terns with 2V = 16 N h = 2,2V = 14 N h = 4 
2V = 12 N h = 6, at J/t = 1.0, h/t = 1.0 and t' = 0. 
presented Ps correspond to singlet pairing in both 
and BCS channel, since at given parameters singlet ] 
ing susceptibilities yield maximum values of Pj> V\ 
Py, is increasing with hole-doping r\ — 2V/,/2V in 1 
channels, the correlation length £ displays different 
havior. In the BCS channel £ first decreases while 
doping decreases from r\ = 2/16 = 0.125 to 4/14 ~ 
and then slightly increases at rj = 6/12 = 0.5. In 
odd channel £ exhibits more than a twofold monot 
increase with doping. 

As was argued in the second paragraph ; 
Eqs. (||,[||l0]), the substantial increase in £ with do 
in the odd channel of the t — J — h model is the cc 
quence of the strongly enhanced spin fluctuations at 
frequencies. This result might indicate that frustra 
due to hole-doping leads to a long-range order in the 
pairing channel. 

In conclusion , we derived composite operators desi 
ing the condensate of superconductor with odd gap J 
generalized t — J model. These composite operator 
ways involve a bound state of Cooper pair and a pr 
combination of magnetization operators. Due to 
presence of an extra spin 1 boson in composite oper 
the parity in odd gap is opposite to that in BCS case, 
odd parity singlet in odd case (see Eq. (|J)) vs. even 
ity singlet in BCS. The BCS and odd pairing correlai 
were investigated numerically for the ID t — t' — J m 
in external staggered field. The enhancement of the . 
correlations coincides with the formation of bound 
pairs and with charge exponent being K p > 1 in 
t — J — h model. It is shown that at small hole doping 
when tf = the strongest odd pairing correlations ar 
ways smaller then strongest BCS correlations. How< 
despite the phase space restrictions, the odd correlat 
increase substantially with hole doping and when t' 
in a frustrated t — t' — J — h model. We argue that do 
and frustration create additional soft spin fluctuat 
what makes odd correlations more pronounced, whe 
BCS correlations are suppressed by doping. 
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FIG. 1. Hole-density correlation function g(r) vs. dist 
r in the t — J — h model, for N = 14, Nh = 4 and two diff< 
values of J/t. 
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FIG. 2. Pairing susceptibilities P vs. distance r in 
t — J — h model. Parameters of the model are the same 
the Fig. (1). Note, that the legend in Fig. (2a) is also 
for Fig. (2b). Pairing susceptibilities are presented onl; 
r > 2 since at smaller distances pairs overlap. 

FIG. 3. The largest pairing susceptibilities in the BCS 
odd channel in the t—t' — J—h model for N = 14, Nh = 4. 
cles present odd singlet pairing, squares present BCS tr 
S z — pairing. 



TABLE I. Values for Pt, and correlation length £ for singlet 
BCS and singlet odd pairing correlations are given for three 
different values of doping rj — 0.125,0.286,0.5 in t — J — h 
model at J/t = 1.0, h/t — 1.0. The pairing correlations in 
triplet channel are smaller. 





BCS 




ODD 




V 


Ps 








0.125 


0.0369 


0.4487 


0.0318 


0.1840 


0.286 


0.0840 


0.2615 


0.0537 


0.2102 


0.500 


0.1924 


0.2864 


0.0838 


0.4100 



